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The photodissociation of the HD molecular ion in intense short-pulsed linearly polarized laser fields is studied using a 
time-dependent wave-packet approach where molecular rotation is fully included. We show that applying a coherent 
superposition of the fundamental radiation with its second harmonic can lead to asymmetries in the fragment angular 
distributions, with significant differences between the hydrogen and deuterium distributions in the long wavelength domain 
where the permanent dipole is most efficient. This effect is used to induce an appreciable isotope separation. 

PACS numbers : 33.80.-b, 33.80.Gj, 33.80.Wz 
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Control of chemical reactivity using two-color laser 
fields has been actively investigated in the past few years. 
The use of a coherent superposition of a fundamental ra- 
diation with one of its harmonics can induce complex in- 
terference effects between different pathways in atomic or 
molecular processes |lj . This method is shown here to be 
a useful tool to control molecular dissociation in the high 
intensity regime. Although several interesting processes 
have been discovered (above threshold dissociation pL[| , 
bond softening J|,|]], vibrational trapping |&-0|), they are 
generally in competition and are not easily disentangled 
in practice. This experimental fact motivated us to deter- 
mine laser parameters which will specifically enhance one 
of these physical processes. Indeed, we have shown that 
in intense laser fields the branching ratios of H^" in the 
different photodissociation channels can be appreciably 
controlled in a (u>, 3u>) experiment Q and that asymmet- 
ric angular distributions of the fragments H + and H(ls) 
can be created in a (lo,2u>) experiment |9. Encouraged 
by the efficiency of this method we decided to explore 
possible asymmetries in the photodissociation of HD + . 

This ion presents an asymmetric electronic distribu- 
tion which is small in the bound molecule but be- 
comes extremely important when the fragments disso- 
ciate. This is because the two lowest electronic states, of 
2 £ + symmetry, asymptotically lead either to H + +D(Is) 
or H(ls)+D + , with an energy difference of 29.8 cm -1 
which is due to the larger reduced mass of the electron 
on the deuterium compared to the hydrogen atom JlO|,[llJ . 
In the ground state, the electron is localized near the deu- 
terium atom while for the first excited state it is on the 
hydrogen. The chemical difference between these two dis- 
sociative channels together with the inherent simplicity 
of this system makes HD + one of the best candidates for 
demonstrating laser-induced control of chemical reactiv- 
ity. 

To simulate the dissociation of the HD + molecular ion 
subjected to an intense and short laser pulse, we solve the 
time-dependent Schrodinger equation using the short- 



time propagator splitting method fjjq| . The molecular 
ion interacts with a classical linearly polarized electric 
field made by a coherent superposition of a fundamental 
radiation of angular frequency ujf and electric field am- 
plitude Ef with its second harmonic of angular frequency 
LJh — 2 ojf and amplitude Eh 

E(t) = /(*) {Ef coB(u f t) + E h cos(w ft i + tp)} e , (I) 

where f(t) — sin 2 (nt/2T p ) is a Gaussian- like pulse shape 
of width T p and total duration 2 T p . The wavefunction 
describing HD + is restricted to include only the two first 
electronic states, with 0c(r, i2) and 4>e(^,R) denoting 
their electronic wavefunctions 



$(R, r, t) = {F G (R, t) G + F E (K, t) <t> E } Xs 



(2) 



R = (i?, </>i{) is the internuclear vector joining the 
hydrogen to the deuterium and r is the electronic coordi- 
nate. The electron spin, described by the Xs function of 
space-fixed projections s = ±1/2, acts here as a specta- 
tor and is uncoupled from the molecular frame. Fg(R, t) 
and Fe (R, t) denote the time-dependent nuclear wave- 
functions associated with each electronic state and con- 
tain the nuclear dynamic information. 

To introduce the angular degrees of freedom of the 
molecule, we expand each nuclear wavefunction in a large 
basis set of spherical harmonics Yn,m{9r! 4>r) (typically 
20 to 40 angular momenta N are necessary to achieve 
convergence 0, depending on the field intensity). We as- 
sume HD + to be initially (t = 0) in a well defined rovibra- 
tional level [v ,N ,M ) of the ground electronic state : 
F G (R,0) = 4> Vo ,N o (R)Y Noi M o (0 R ,(f>R) and F E (H,0) = 0. 
For a linearly polarized electric field the initial M quan- 
tum number is conserved, and the wavefunctions Fq (R,t) 
and Fe (R,t) are expanded on a set of spherical harmon- 
ics with fixed azimuthal quantum number 

F G{E) (R,t)=J2 F G(E),N(R,t)Y NMo (0R^R) ■ (3) 

N 
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In contrast to the H^" case 0j, this expansion in the 
spherical harmonics basis includes both parities for each 
electronic state, due to the permanent dipole moment in 
HD + which allows transitions within a given electronic 
state. Indeed, the radiative interaction directly couples 
the nuclear wavefunctions of components N <-> N ± 1. 
The corresponding coupling matrix elements, evaluated 
in the length gauge, take the same form as for H^~ 

VZ%_ l{ R, t ) = ^(R) J {2N N _\- { ^ + 1} E(t) , (4) 

except for the expression of the different dipole moments 
fJ-a/3 (R) between a pair of electronic states denoted by the 
subscripts a and (3= G or E. To take the non-degeneracy 
effects into account, and thus present the correct disso- 
ciation behavior, we use the transformed Hamiltonian 
defined by Moss and Sadler , and employ the coupled 
states approach of Carrington et al [ pj| . We expand the 
resultant ground G and first excited 4>e electronic wave- 
functions in terms of the pure lscr ff and 2po~ u solutions of 
the Born-Oppcnheimcr Hamiltonian 

4> G = a{R) <pl sag +b(R) (f> 2 pa u r~oo +01s( D ) (z\ 

$ E = -b(R) 4>u„ 9 +a(R) fop« u R ^-MH) { > 

where a ~ 1 and b ~ for short internuclear distances. 
At large distances, both a and b approach 1/V2 and the 
(j)G and <pE wavefunctions correlate with the Is atomic 
orbitals centered on the deuteron or on the proton re- 
spectively. On the way to dissociation, the electron 
can no more "jump" from one fragment to the other, 
and the electronic dipole moment /ieg(R) vanishes at 
long distance. But in this limit, the internal couplings 
Hgg(R) and hee(R) diverge as —2R/3 and R/3 respec- 
tively The non-adiabatic couplings induced by 
the non-diagonal elements of the d 2 /dR 2 kinetic energy 
operator in the (G,E) electronic basis are also included 
in the time propagation using an additional basis trans- 
formation. 

At the end of the propagation, the partial nuclear 
wavefunctions FG(E),N(R,t — 2T p ) are projected on the 
corresponding field-free continuum nuclear wavefunctions 
to get the dissociation probabilities. One thus obtains the 
energetic and angular distributions of the fragments , 
to be compared with experimental data when available. 

A simplified four-state model can be derived from the 
above multistate treatment in analogy with the two-state 
model commonly used for the Hj ion in intense laser 
fields if one assumes that the molecular ion is totally 
aligned by the laser electric field along its polarization 
axis. This alignment has been demonstrated experimen- 
tally by different groups jl3| and results from the optical 
pumping of the rotational quantum number N S. To 
take the effect of the permanent dipole moments /x GG (i?) 
and [1ee{R) into account, we associate two nuclear wave- 
packets with each electronic potential, representing the 



collective sum of the odd (o) and the even (e) partial 
wave components in Eq.(3). The total wavefunction in 
Eq.(2) can then be expressed simply as follows 

$ = {{F G ,o + F G>e ) <t> G + (F E .o + F E . e ) <Mx s • (6) 

We propagate an initial (t — 0) wavepacket taken as a 
pure w vibrational state of the ground electronic po- 
tential. The /x GG (i?), (J-ee{R) and [Ieg(R) dipole mo- 
ments are included in the potential matrix as off-diagonal 
terms that couple nuclear wavefunctions of opposite par- 
ity. This simple model saves a lot of time and memory 
in computation and has been checked to give reliable re- 
sults. The angular distribution of the H + fragment can 
be obtained from the projections C Gi0 (e) and C Gje (e) of 
the odd and even G nuclear wavepackets at the end of the 
pulse on the G field-free continuum nuclear wavefunction 
at the energy e by 

P H+ M = 0°) = i|C G . (e)-C G , e ( £ )| 2 
P H +M = 180°) = ±|C G , (e) + C G , e ( £ )| 2 , 1 1 

where 8 = 0° and 9 = 180° denote the forward and 
backward ion signals pjj^]. The directional dissociation 
probabilities are obtained by energy integration 

Ph+ (0 = 0° or 180°) = J P H +(e, = 0° or 180°) de , (8) 

and the total dissociation probability Pjj+ m t he ground 
electronic state is just the sum of these two directional 
probabilities. Similar equations (with opposite signs in 
Eq.(7)) relative to the excited part E of the wavepacket 
give the partial and total productions of D + ions. 

The most remarkable effect on the angular distribu- 
tions of the HD + photofragments is obtained using a 
coherent superposition of the second (or any even) har- 
monic with a fundamental laser frequency. In this case a 
given energy of dissociation e can be reached by absorb- 
ing odd as well as even number of photons. Thus, for 
example, Cg,o{£) and C G! e(e) can simultaneously con- 
tribute to the probabilities in Eq.(7), and the cross- term 
Re{C G Q (e) C Gi e(e)} induces different dissociation prob- 
abilities in the 6 = 0° and 6 = 180° directions. This 
effect is similar to the asymmetric angular distributions 
of photoelectrons observed in two-color high intensity 
atomic photoionization El] , or to the asymmetric proton 
angular distributions in photodissociation How- 
ever, what is unique here is that we can also induce dif- 
ferent asymmetries in the distributions of the H + and 
D + photofragments and control the spatial separation 
of these isotopes. This effect requires long wavelength 
lasers, such as the CO2 laser case shown in figure 1 us- 
ing Xf = 10.6/xm (TiiOf = 970cm _1 ) and Ah = 5.3/^m 
(hw h = 1940 cm- 1 ). 

This figure presents the fully converged 3D results ob- 
tained from an expansion using up to 40 coupled par- 
tial waves in each electronic potential (see Eq.(3)), for 
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an initial (v Q — 0, N a — 0, M a = 0) rovibrational state 
and using peak intensities // = 5 x 10 13 W.cm -2 and 
Ih = 2 x 10 13 W.cm~ 2 to evaluate the electric fields in 
Eq.(l). The pulse duration is 225 fs, corresponding to 13 
optical cycles of the second harmonic. For each chosen 
phase (p the forward and backward dissociation proba- 
bilities for a particular fragment are obtained by angular 
integration of the 3D distributions over 0° < 9 < 90° and 
90° < 9 < 180° respectively where 9 denotes the angle 
between the dissociation and the polarization axis. 



and N = 0, 1, 2 and M = 0, ±1. This coherent control 
should thus persist even after averaging over an initial 
distribution of the HD + states, as expected in a real ex- 
periment. 



FIG. 1. Dissociation probabilities of HD + in the two first 
electronic states leading to the H + + D(ls) (dashed line) and 
H(ls) + D + (solid line) fragments as a function of the phase 
shift ip between the two laser harmonics. The integrated 
forward (9 ~ 0°) probabilities are shown in (a), and the 
integrated backward (9 w 180°) probabilities are shown in 
(b), for a laser pulse of 225 fs duration made by the coher- 
ent superposition of the two wavelengths A/ = 10.6 /im and 
\h = 5.3 /im with the intensities If — 5 x 10 13 W.cm -2 and 
I h = 2x 10 13 W.cm' 2 . The HD + molecular ion is initially in 
the (v = 0, N = 0, Mo = 0) rovibrational state. 

The polar plot in figure 2 for the specific phase ip = 
confirms that the angular distributions are closely peaked 
along the polarization axis e. For this phase shift ip = 0, 
essentially 100% of the dissociated ions emitted in the 
forward direction are expected to be H + , and the D + 
ions are almost only ejected in the backward direction. 
Obviously this gives a very efficient mechanism for spa- 
tially separating the isotopes. 

In this dissociation scheme, a neutral atom H(ls) or 
D(ls) always appears in the opposite direction to the 
other ion (D + or H + ). A significant ionization of the 
neutral atoms on the way to the detector should thus 
just enhance the ion signal, and increase the efficiency of 
the separation process. If the phase shift is changed to 
<p = 7r the direction of ejection of the H + and D + ions 
are simply reversed, due to a change in the sign of the 
cross-term in Eq.(7). It is just near tp = tt/2 that the 
isotopes are not separated. 

The explorations we have made show that this sep- 
aration effect is robust enough to survive experimental 
uncertainties in the intensities of both harmonics and in 
their phases. Increasing the pulse length from 225 fs to 
675 fs has a negligible effect on the angular asymmetry 
for </? = 0. We further tested the experimental feasibil- 
ity by performing the same calculations for a wide range 
of initial vibrational and rotational states and found the 
same qualitative behavior as in figure 1 for v = 0, 1, 2, 3 



FIG. 2. Polar plot showing the angular distributions of the 
H + and D + ions with the same conditions as in figure 1 for 
the specific phase tp — 0. Note the better alignment of D + , 
due to a more efficient rotational pumping |)| in the excited 
electronic state. 

The most sensitive parameter for achieving coherent 
control of the isotope separation in HD + photodissocia- 
tion is the laser frequency LOf in Eq.(l). Indeed, twice 
the CO2 laser photon energy (1940 cm -1 ) is compara- 
ble to the vibrational spacing in HD + ground electronic 
state. A choice of laser frequency in this range greatly 
enhances the efficiency of the ^igg(R) &ud fJ-EE(R) per- 
manent transition dipoles. The variations of the forward 
and backward dissociation probabilities as a function of 
the fundamental wavelength used in this two-color study 
are shown in figure 3 for the most favorable phase tp = 0. 
These ID calculations are presented for the initial level 
v = 3 using the same intensities as in figure 1. The re- 
markable separability seen in figure 2 persists over the 
full range of frequencies 500 cm -1 <huif< 2500 cm' 1 . 

It is interesting to note that there are two distinct 
regions of enhanced probability of photodissociation in 
figure 3 which we interpret as low and high frequency 
mechanisms. The first one is analogous to the static 
field tunneling regime found in strong field photoioniza- 
tion and photodissociation , while the second is closer 
to the Floquet multiphoton regime. The long wavelength 
(low frequency) mechanism, where appreciable dissocia- 
tion occurs during a single optical cycle |15|, favors the 
spatial separation of the isotopes. In this regime, path- 
ways induced by the permanent dipoles /j,gg and Hee, 
which have opposite signs, manifest interferences with 
those due to the heg transition dipole, and it allows us 
to achieve this remarkable control. On the contrary, the 
fragmentation dynamics in the short wavelength (high 
frequency) domain is characterized by two features which 
destroy this interference effect. Firstly, because nothing 
significant happens during one optical period, the dissoci- 
ation dynamics reflects behavior "averaged" over the op- 
tical cycle. Secondly, unfavorable Franck-Condon factors 
between nuclear wavefunctions differing by the energy of 
one photon within a given electronic state reduce the ef- 
fect of permanent dipoles. In fact, since the (J-eg coupling 
dominates the dissociation in this high frequency limit, 
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the dynamics of HD + is very similar to the one of and 
no appreciable difference between the H + and D + partial 
probabilities can be induced in this regime, as observed 
by Sheehy, Walker and DiMauro [ fl6| |. 



FIG. 3. ID-dissociation probabilities of HD + (v = 3) in 
the forward (a) and backward (b) directions leading to the 
H+ + D(ls) (dashed line) and H(ls) + D+ (solid line) frag- 
ments as a function of the fundamental wavelength. The 
phase of the second harmonic is fixed to <p — 0. The pulse 
duration and laser intensities are the same as in figure 1. 

In this paper we have predicted a new and efficient 
mechanism for separating different isotopic products in a 
fragmentation process. Unlike previous examples of co- 
herent control which only controlled the distribution of 
all the charged particles together, we are now able to dis- 
tinguish the fragments of different masses. This scenario 
should extend far beyond the specific example of HD + 
and can be applied to separating products of a different 
chemical nature as well. Tailoring quantum interferences 
induced by two-color coherent light opens a promising 
avenue for controlling the dynamics of quasi-symmetric 
systems, and applications of primary importance can be 
expected in enantio-differentiation for instance fTij ]. 
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